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ABSTRACT 



Context. Although water is an essential and widespread molecule in star-forming regions, its chemical formation pathways are still 
not very well constrained. Observing the level of deuterium fractionation of OH, a radical involved in the water chemical network, is 
a promising way to infer its chemical origin. 

Aims. We aim at understanding the formation mechanisms of water by investigating the origin of its deuterium fractionation. This 
can be achieved by observing the abundance of OD towards the low-mass protostar IRAS 16293-2422, where the HDO distribution 
is already known. 

Methods. Using the GREAT receiver on board SOFIA, we observed the ground-state OD transition at 1391.5 GHz towards the low- 
mass protostar IRAS 16293-2422. We also present the detection of the HDO ln-Ooo line using the APEX telescope. We compare the 
OD/HDO abundance ratio inferred from these observations with the predictions of chemical models. 

Results. The OD line is detected in absorption towards the source continuum. This is the first detection of OD outside the solar 
system. The SOFIA observation, coupled to the observation of the HDO ln-Ooo line, provides an estimate of the abundance ratio 
OD/HDO ~ 17-90 in the gas where the absorption takes place. This value is fairly high compared with model predictions. This may 
be reconciled if reprocessing in the gas by means of the dissociative recombination of H 2 DO + further fractionates OH with respect to 
water. 

Conclusions. The present observation demonstrates the capability of the SOFIA/GREAT instrument to detect the ground transition 
of OD towards star-forming regions in a frequency range that was not accessible before. Dissociative recombination of H2DO + may 
play an important role in setting a high OD abundance. Measuring the branching ratios of this reaction in the laboratory will be of 
great value for chemical models. 

Key words. OD - SOFIA - GREAT instrument - IRAS 16293-2422 - deuterium fractionation - astrochemistry 
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1. Introduction 

Water is an essential molecule in star-forming regions because 
it is in certain co nditions one of the main gas coolants (e.g. 
iNisini et all l2010h . and it is the most abundan t constituent o f 
the icy mantles frozen onto dust grains (e.g. Gibb et al. 2004). 
Despite this fame, water is still mostly a mystery for astro- 
chemists: its formation pathways are not very well constrained. 
Water can form in the gas phase at high temperatures, or on the 
dust grains via surface chemistry. Surface chemistry in particu- 
lar is believed to play an important role in the cold and dense 
clouds where stars form. But the details of these processes are 
mostly unknown. The reaction network involves hydrogenation 
of oxygen atoms (O) and O3 to form OH, which can in turn react 
with H or H2 to form H2O, as well as succes sive hydrogena- 
tions of O?, leading to HOOH and water Ce.g. lTielens & Hagenl 
1982). This reaction network is at the moment t he subject of 
intense laboratory studies (e.g., lOba et all 120091: lloppolo et alJ 
2010; ICuppen et alJ l2010h . Some new observational insights 
into the water formation pathways have been achieved with the 
first detection in spa ce of hydrogen peroxi de (HOOH) towards 
the p Oph A cloud (tBergman et al.N201lh . HOOH is the pre- 
cursor of water in the hydrogenation of O2 on grain surfaces. 
Detailed coupled gas-grain chemical modelling of this detection 
was able to simultaneously reproduce the observed abundances 
of for maldehyde, met hanol, O2 and hydrogen peroxide in this 
source (iDu etal.l2012l) . validating the importance of grain chem- 



istry processes in the formation of HOOH and as a consequence 
in water formation. 

We propose here yet another, complementary, approach for 
the study of water formation. The level of deuterium fractiona- 
tiorQ in molecules has been shown to be a good signpost of their 
formation conditions. Molecules forming in cold and dense en- 
vironments are observed to show a high level of fractionation, 
which even led to the detection of tr iply-deuterate d molecules 
such as ND 3 (Ivan der Tak et al.ll2002l) and CD 3 OH dParise et all 
2004). The observation of the deuterium fractionation in a 
molecule and its molecular precursors has proven to be very 
useful in constraining the relative importance of formation path- 
ways. For example, the observation of formaldehyde (H2CO) 
and isotopologues (HDCO and D2CO) as well as methanol 
(CH3OH) and isotopologues (CH2DOH, etc.) in a sample of 
low-mass protostars has shown that methanol is mainly formed 
on dust surfaces by hydrogenation of CO, whereas a signifi- 
cant fraction of fo rmaldehyde must be formed in the gas phase 
(Paris e et alj|20 06). Another promising example is the study of 
the NH and ND radicals with the aim to understand the forma- 
tion of NH3 and isotopologues (Bacm ann et al.ll2010h . 

Water formation on dust surfaces is believed to proceed via 
the hydrogenation of either OH or HOOH, the latter case lead- 
ing to the production of an additional OH. The abundance of 



1 Enrichment in deuter ium in a mol ecule, compared with the cosmic 
D/H ratio (~ 1.5 x 10~ 5 , Linsky 2003h . 
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Fig. 1. Lowest energy levels of the OD radical. The arrow shows 
the observed line at 1391.5 GHz. Each rotational level on this 
figure is split to first order into a A-type doublet (/ = +1), each I 
level is then split again into hyperfine sublevels (see Fig. [5]). 

H2O and OH should be therefore tightly linked if grain processes 
are dominant. The deuterium fractionation of water in hot cori- 
nos (i.e. in the gas enriched from evaporation of dust mantles 
in low-mass protostars) is found to be anomalously lo w in com- 
pariso n with other molecules (e.g., Parise et al. 2005; Liu et al. 
|2011|) . and the reason for this is not yet understood. Therefore, 
observations of the OD/OH ratio in these same regions will 
help to understand the peculiarity of water chemistry in young 
stellar objects. In the present study, we target OD towards the 
low-mass protostar IRAS 16293 -2422 (hereafter IRAS 16293), 
where t he distribution of d euterated water has been previously 
studied dParise et al.ll2005l) . We note that the OD/OH ratio is 
predicted to be relatively high if it is formed exclusively in 
the gas phase (mor e than -0.1 at T=10K and «# = 10 4 cm" 3 , 
iRoberts et ai]|2002l) . 

OD has been searched to date unsuccessfully towards the 
Galactic Centre, by means of the observation of the AF = 0A- 
doublet transitions of the 2 n3/2, J=3/2 state around 310 MHz 
(lAllen et al.lll974l) . This non-detection may now be explained 
by the lower D/H ratios observed in general towards the Galactic 
Centre, due to astration (destruction of deuteron nuclei in stellar 
interiors). OD has been recently detected in the comet C/2002 
T7 (LINEAR) via the co-addition of 30 lines of its ultraviolet flu- 
orescence spectrum with a ratio [OD]/[OH] = (3.5+1.0) x 10~ 4 
(Hutseme kers et al.l2008l) . These observations, interpreted as the 
result of photodissociation of HDO and H2O, lead to an estimate 
of the water fractionation of the same order as in other comets, 
a factor of 10 lower than observed in hot corinos. 

The present paper is organised as follows. Section |2]briefly 
describes the spectroscopy of OD. Section [3]presents the obser- 
vations, which are analysed in Sect. [4] Results are discussed in 
Sect. 



2. OD spectroscopy 

Figure Q] shows the lower energy levels of the OD molecule. 
Owing to the unbalanced electronic spin and orbital angular mo- 
menta, the spectrum of the OD radical shows A-type doubling, 
similar to that of the OH radical. Each rotational level is split 
into a A-type doublet. Each component of the doublet is fur- 
ther split by the magnetic hyperfine structure (hfs) due to the 
deuterium atom into components characterised by the total an- 
gular momentum F (=J+I). F takes values of J-l, J and J+l 
(Dousmanis et al. 1955). The selection rule of AF=0,+1 leads to 
six hfs lines (Fig. |5) . We use here t he transitions tabulated in the 
JPL database^ (Pic kett et al.| [l998). These predictions are based 



on the laboratory measureme nts of iBeaudet & Povnterl Q1978) 
and lBrown & Schuber3 (ll982l) . 



3. Observations 



Using the GREAT instrument (H evminck & al]l2012l) on board 
SOFIA, we observed the OD 2 n 3/2 J=5/2^3/2, 1=-1^+1 
line multiplet towards the low-mass protostar IRAS 16293 at the 
coordinates a(2000)=16 h 32 m 22 s :90, <5(2000)=-24°28'36'.'3. The 
double-sideband (DSB) receiver was tuned in the upper sideband 
at the frequency 1391494.71 MHz of the strongest hfs compo- 
nent. The receiver was connected to a digital FFT spectrometer 
dKlein et alJ2012h providing a spectral resolution of 0.04 km s" 1 , 
and a bandwidth of 1.5 GHz. Two setups were observed with 
different LSR velocities (-2kms _1 and +8kms~') to disentan- 
gle any contamination from the image sideband. We used the 
double-beam switch mode chopping in E-W direction with an 
amplitude of 70" and a phase time of 0.5 sec. The atmospheric 
transmission was fit with special care during the calibration pro- 
cess to remove a residual telluric line at 30 km/s. Part of the data 
suffer from a pointing error in the second phase of our double- 
beam observations. Because the OD line is observed in absorp- 
tion against the continuum of the source, which is likely compact 
at this frequency compared with the SOFIA beam (20'.'4 at the 
OD frequency), a pointing offset should only result in a differ- 
ent coupling of the signal with the telescope beam. We therefore 
scaled all ON-OFF cycles to the brightest continuum level and 
averaged them. The resulting spectrum is obtained with 1 1 min 
on-source integration time. An additional standing-wave in the 
signal was removed by killing the appropriate peak at ~0.05 
MHz" 1 in the Fourier plane of the signal. The resulting spec- 
trum is presented in Fig. [2] The beam efficiency is taken to be 
0.54 dHevminck & alJl2012h . The OD feature clearly appears in 
absorption, with a peak signal-to-noise ratio of 6 at the resolution 
of 0.79 km s" 1 , the rest of the band being free of any emission or 
absorption. The DSB-continuum is 2.7+0.2 K at the same reso- 
lution. 




-50 



Velocity (km/s) 



Fig. 2. OD J=5/2— >3/2, 1=-1— >+l transition observed towards 
IRAS 16293. The reference frequency is that of the main hfs 
component. The spectrum is smoothed to 0.79 km s" 1 resolution. 
The continuum level is twice the true continuum because of DSB 
calibration. 



We also present here observations of the HDO 1 i,i-Oo,o line 
obtained with the CHAMP + receiver at APEX. The observa- 
tions were performed on 2010, September 6 and 7, under very 



http://spec.jpl.nasa.gov/ 



GREAT is a development by the MPI fur Radioastronomie and 
the KOSMA/Universitat zu Koln, in cooperation with the MPI fur 
Sonnensystemforschung and the DLR Institut fur Planetenforschung. 
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Fig. 3. HD01i,i-Oo,o and OD line profiles observed towards 
IRAS 16293. The HDO spectrum was shifted by -2.9 K along the 
y-axis. The OD spectrum was smoothed to 0.39 km s -1 . The hy- 
perfine intensity structure of the OD line is shown as red bars. 
Note that the OD observation is DSB while the HDO observa- 
tion is SSB. The green curve shows the two-Gaussian fit of the 
HDO profile (see Table[fl. 

stable atmospheric conditions (PWV = 0.35 -0.40 mm). The 
CHAMP+ high-frequency channel was tuned to 893.639 GHz, 
while the low-frequency channel was tuned on CO(6-5). The fo- 
cus was checked on Venus, and the local pointing was checked 
every 1 to 1.5 hours with CO(6-5) line pointing on NGC6334I. 
Pointing corrections were found to be of the order of 3 " . The an- 
gular resolution of the HDO observations is 7 " . The receiver was 
connected to the AFFTS, providing a resolution of 0.06 km s -1 , 
and a bandwidth of 1.5 GHz. The continuum is 2.9+0.2 K (T m & 
scale) at resolution 0.37 km s -1 . 



HDO line has a similar value in the 50- 100 K range. We can 
therefore reasonably assume that both absorptions occur in the 
same gas, and that the excitation temperature of the lines will 
be extremely low, as the density of the IRAS 16293 envelope is 
several orders of magnitude lower than this critical value. 

We assume that each hyperfine component of OD has 
an intrinsic line width similar to that of the HDO absorp- 
tion (lkms -1 ), and we model the global absorption using the 
XCLASS softwareQ. The OD line is close to saturation, as al- 
ready visible in Fig. [3] where the depth of the absorption is about 
half that of the DSB continuum (i.e. the line absorbs the total- 
ity of the continuum in the signal sideband, leaving the image 
sideband continuum unabsorbed). The saturation remains lim- 
ited (t < 4), as no line broadening (except for that due to hyper- 
fine structure) is observed. The signal-to-noise ratio of the ob- 
servation is not high enough, however, to infer the optical depth 
from the comparison of the depth of the three main hyperfine 
components. The OD and HDO column densities inferred from 
the XCLASS modelling averaged over the extent of their respec- 
tive continuum for different assumed T ex are tabulated in Table|2] 
The uncertainties are estimated by comparing by eye the obser- 
vation with the synthetic model. Because both lines are seen 
in absorption against the source continuum, and the absorbing 
layer is likely to be more extended than the continuum in each 
case, we can compare the two column densities without beam 
dilution hindrance. With these assumptions, the abundance ratio 
OD/HDO is found to be in the range 17-90, if we assume that 
the OD and HDO absorptions come from the same region of the 
envelope, and are described by the same excitation temperature, 
uniform along the line-of-sight. 



4. Analysis 

Figure [3]presents the OD spectrum observed with SOFIA (black 
curve), as well as the HDO line observed towards the same 
source with APEX. The OD spectrum is characterised by a broad 
absorption (2.9 km s -1 , see Table [TJ, compared with the HDO 
line (l.Okms -1 ). This broadening is due to the hyperfine struc- 
ture of the OD line. An hfs fit using the CLASS softwareQ pro- 
vides an intrinsic linewidth of 1.3+0.6 kms , which is consis- 
tent with the width of the HDO absorption within the error bars. 



Table 1. Fitting results. 



Molecular line 


Fit type 


/T,„z,dv 
(Kkms" 1 ) 


FWHM 

(kms" 1 ) 


Vfc, 

(kms" 1 ) 


HD0 1 u -0 ,o 


two- Gauss 


+4.4+0.6 
-2.4+0.3 


5.9+1.0 
1.0+0.1 


4.5+0.3 
4.2+0.1 


OD 5/2->3A -1-H-l 


one-Gauss 


-4.1+0.5 


2.9+0.3 


- 


OD 5/2->3/2, -1-M-l 


hfs 


- 


1.3+0.6 


4.2+0.2 



Collision rate coefficients for the OD molecule are not avail- 
able. However, an estimate of the critical density of the upper 
level of the transition can be obtained by looking at the critical 
density of the corresponding line of OH (at 2514 GHz). Using 
the values of the Einstein coefficient for spontaneous emission 
and the collision rates tabulated in the LAMDA database^, we 
find that the critical density of the corresponding OH upper level 
is abound 10 10 cm 3 at 60 K. The critical density of the 894 GHz 



4 http://iram.fr/IRAMFR/GILDAS/ 

5 http://www.strw.leidenuniv.nl/~moldata/ 



Table 2. 
elling. 



Column densities inferred from the XCLASS mod- 



T„ (K) N OD (cm- 2 ) 



N HD0 (cm- 2 ) OD/HDO 



2.7 
5.0 
10.0 



(3.5+1.5) xlO 13 
(5.0+2.0) xlO 13 
(1.0+0.3) xlO 14 



(6.0+1.5) xlO" 60: 
(1.2+0.3) xlO 12 45: 
(4.0+1.0) xlO 12 27: 



30 
20 
10 



5. Discussion 

Ideally, we would like to measure the OD/OH ratio in the en- 
velope of IRAS 16293 to compare it on the one hand with the 
HDO/H2O ratio in the same environment, with the aim to con- 
strain astrochemical models describing the synthesis of water 
and its deuterium fractionation, and on the other hand with the 
OD/OH ratio in more evolved environments, e.g. in comets, to 
infer the amount of chemical reprocessing during the formation 
of a solar-type star. To our knowledge, no column density for 
OH in the absorbing envelope of IRAS 16293 has been reported 
so far. Although Herschel/HIFI could provide a high-resolution 
observation of the OH 2 n 1/2 J=3/2 -> 1/2 line at 1.83 THz, this 
line has an upper energy level too high to be directly compared 
to the OD line presented here. A more promising perspective in 
the future will be to observe the corresponding OH line ( 2 H}/2 
J=5/2 -> 3/2) at 2.5 THz with the upgraded GREAT instrument. 
In the meantime, we discuss here the value of the OD/HDO 
column density ratio inferred in Sect. |4] In the absence of 
a differential fractionation in OH and water, one expects the 
OD/HDO ratio to be close to the OH/H2O abundance ratio. 



http://www.astro.uni-koeln.de/projects/schilke/XCLASS 
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Fig. 4. Left panel: Predictions of the OH/H2O ratio from the iDu et all (120 1 2|) standard model, for constant (T, hh) conditions. Right 
panel: Models with con stant and time-d ependent temperature. "Low Eb" models use the binding energies for species on the grains 
of the standard model of lDu et all J2012I) while "high Ef," models use the higher binding energies for species with an OH bond from 
Garrod & Herbst (2006). 



lDu et alJ d2012h have presented a chemical model that in- 
cludes processes in the gas-phase and on the grains, to describe 
the formation of water. We use the same model here, consider- 
ing constant physical conditions, in the range T — 10-30 K and 
«h - 10 4 - 10 6 cm -3 (see Fig. |4] left panel). Except at the very 
early times (t< 10 yr), the gas-phase OH/H2O ratio is always 
predicted to be fairly low compared with our observed OD/HDO 
ratio. The closest agreement happens for times later than 10 4 - 
10 5 yr, when the OH/H2O can reach a maximum value of 5.7. 

The models stated above assume constant temperature and 
density conditions. However, during the formation of a star, the 
envelope is heated by the central protostar, and a model where 
the temperature evolves with time may be more appropriate. 
Figure|4](right panel) shows some results with a t ime-dependent 
temperature based on the model of lDu e t al. (2012). For the mod- 
els with protostellar heating, T increases between 12 K and 50 K 
linearly with time, either in the 1-3 10 5 yr ("fast" heating) or in 
the 1-10 10 5 yr range ("slow" heating). We also test different sets 
for the desorption energies of species from the grains. None of 
the models seem to be able to explain the high inferred OD/HDO 
ratio. There could be two reasons for this discrepancy: either the 
radiative transfer we used is too simplistic (in which case OD 
collisional rates are highly needed to go beyond this treatment), 
or the OH radical gets much more fractionated than water. This 
latter case ma y arise because the OH + D exchange reaction 
is exothermic (Roberts & Millar 2000). The relative fractiona- 
tion of OH with respect to H2O could be even more enhanced 
if the gas-phase dissociative recombinations (DR) of H30 + and 
H2DO + have non-trivial branching ratios for the formation of 
water and hydroxyl. Our model based on the UMIST06 rate as- 
sumes a produced OH:H20 ratio of 3: 1, consistent with the work 
of lJensen e t al. (2000), who studied the DR of all deuterated iso- 
topologues of H30 + but H2DO + . Their results concerning the 
DR of HD 2 + shows that both HDO and D 2 are produced in 
similar amounts, whereas the OD:OH ratio is around 3.4, which 
may explain our high observed OD/HDO ratio. Experiments on 
the DR of H2DO + are needed, however, before drawing a firm 
conclusion. 



6. Conclusions 

This Letter presents the first detection of OD outside the solar 
system, and demonstrates the capability of the SOFIA/GREAT 



instrument for breakthrough observations of spectral lines at 
high spectral resolution in a frequency band previously not 
covered (the OD line presented here falls into the frequency 
gap of the HIFI receiver on board Herschel). The high in- 
ferred OD/HDO ratio in the envelope may be the signpost of 
gas-phase reprocessing through dissociative recombination of 
H 2 DO + . After the commissioning of the GREAT 2.5 THz (M- 
)channel, parallel observations of the OD and the OH ground- 
state transitions could be performed routinely towards, a.o., star- 
forming regions. The observation of these radicals is expected to 
provide valuable constraints on the formation and fractionation 
of water. 
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Fig. 5. Hyperfine structure of the 1391.5 GHz transition. The spacing between the different hyperfine levels was increased for 
visibility, and the overall separation between the J=5/2 and 3/2 levels is not to scale. The transitions J=5/2— »3/2, 1=+1— »— 1 at 
1390.6 GHz, difficult to observe from SOFIA because of residual atmospheric absorption, are not shown. 



